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The oxide pyrochlore Bi2Ti20eO' is known to be associated with large displacements of Bi and 
O' atoms from their ideal crystallographic positions. Neutron total scattering, analyzed in both 
reciprocal and real space, is employed here to understand the nature of these displacements. Rietveld 
analysis and maximum entropy methods are used to produce an average picture of the structural 
non-ideality. Local structure is modeled via large-box reverse Monte Carlo simulations constrained 
simultaneously by the Bragg profile and real-space pair distribution function. Direct visualization 
and statistical analyses of these models show the precise nature of the static Bi and O' displacements. 
Correlations between neighboring Bi displacements are analyzed using coordinates from the large- 
box simulations. The framework of continuous symmetry measures has been applied to distributions 
of 0'Bi4 tetrahedra to examine deviations from ideality. Bi displacements from ideal positions 
appear correlated over local length scales. The results are consistent with the idea that these 
nonmagnetic lone-pair containing pyrochlore compounds can be regarded as highly structurally 
frustrated systems. 



PACS numbers: 61.05.&n, 61.43. Bn, 75.10.Jm 
I. INTRODUCTION 

Magnetic oxides with the A2-B2O7 pyrochlore struc- 
ture have been the subject of intense study. The lattice 
of corner-connected tetrahedra of A atoms hinders coop- 
erative magnetic ordering, and when the A atom spins 
are Ising, an ice-like ground state is produced, for exam- 
ple in Dy 2 Ti 2 7 and HoiTiaOT™ Other recent devel- 
opments in oxide pyrochlores include superconductivity 
in osmium compounds,^ the formation of polar metallic 
states with unusual phonon modest, and the suggestion 
of chiral magnetic ground statesP The analogy between 
Ising spins and vector displacements of cations within 
their coordination polyhedra has led to the suggestion 
that polar ordering may be similarly frustrated on the 
pyrochlore latticeP Thus in pyrochlore Bi2Ti20gO', it 
is known that the Bi 3+ atoms, usually predisposed to 
off-centering within their coordination polyhedra, dis- 
play incoherent displacements permitting the average 
structure to remain cubicP This is in sharp contrast 
to the Bi 3+ -containing perovskites BiMnOa and BiFe03 
where the lone-pair active A-site produces polar, non- 
cubic ground states.^ Some signatures of these inco- 
herent displacements are seen in measurements of heat 
capacity of Bi2Ti20gO' and related compounds at low 



temperatures.^ 

Significant advances in describing frustrated, ice-like 
behavior in magnetic pyrochlores have been made when 
atomistic models arc utilized to describe the interactions 
between individual A sites. For example, atomistic sim- 
ulation of magnetic spin ordering in Dy2Ti207 leads to a 
picture of localized, un com pensated spins connected by 
strings of ordered spinsP^ This picture is local: spin be- 
havior is driven by connectivity, geometry, and pairwise 
exchange interactions. However, it also agrees with bulk 
thermodynamic measurements: the heat capacity of the 
ensemble average agrees with experimental observations. 
We show here that experimental modeling over multi- 
ple length scales, both atomistic and averaged, affords 
a view of the pyrochlore Bi2Ti2 0eO' where Bi displace- 
ments, rather than spins, form a frustrated network on 
the A sublattice. 

Bi2Ti20gO' [structure in Fig.fija)] is written thus to 
emphasize the two sublattices: one of corner-sharing 
TiC>6 octahedra, and the other of corner-sharing 0'Bi4 
tetrahedra. While the TiOg sublattice is rigid in mod- 
els of the average structure, Bi atoms are suggested 
to displace 0.4 A normal to the linear O'-Bi-O'bond 
in an uncorrelated mannerP^ First-principles calcula- 
tions on Bi2Ti20eO' predict Bi displacements but these 
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are perforce associated with non-cubic symmetries.^ 
Diffuse intensity in electron diffraction patterns of re- 
lated compounds including Bi 2 Ru 2 07, Bi 2 lnNb0 7 , and 
Bi 2 ScNb07 may indicate short-range correlations in the 
Bi displacementsP^HEl If Bi displacements cooperatively 
order with each other, they must do so only over short 
ranges. Crystallographic analysis based on Bragg scat- 
tering leaves a void in the ability to probe such short- 
range order, as analyses are predicated on the existence 
of long-range order. Consequently, studies of displacive 
disorder on the A site via Rietveld refinement or Fourier 
maps can produce a model of the average electron or nu- 
clear distributions, but each A site has an identical cloud 
of intensity.^ 

We investigate models where the correlated motion of 
atoms on the A sites reproduces the atomistic, pairwise 
distances between individual atoms. This description is 
provided by an appropriate Fourier transform of the total 
scattering function S(Q) to pro vide a normalized pair dis- 
tribution function (PDF) In this study, the PDF and 
the Bragg profile are used as experimental constraints in 
a large-box (11,000 atom) model of Bi 2 Ti 2 6 0' to ob- 
tain, using reverse Monte-Carlo (RMC) analysis, a con- 
sistent picture of the the coordination tendencies of all 
atoms. Many of these models are combined and used as 
a set of atomic positions for further analysis. RMC com- 
pares the experimental and computed (based on atom po- 
sitions in the simulation box) D(r) and S(Q) while ran- 
domly relaxing atomic positions. The method is similar 
to Metropolis Monte Carlo, except that the fit to data y 2 , 
instead of a potential energy function, is minimizedP^HSB 

There are many approaches to describing the behav- 
ior and correlations of atomic positions as obtained from 
large-box models of structure. Some examples include 
the use of quadrupolar moments of octahedra to de- 
scribe LaMnC>3p2l of contour plots of bond angles in 
cristobalitep^ and the use of bond valence analysis to 
obtain valence st ates f rom a statistical analysis of metal- 
oxygen positions! 24 * 25 ^ Here we analyze the local geom- 
etry using simple metrics, then present the continuous 
symmetry measures (CSMJpHHor polyhedra from RMC 
simulations. The CSM model provides a quantitative 
measure of a polyhedron's distortion, in the form of a 
"distance" from ideality. A key advantage of CSM is its 
ability to compare shapes in different compounds.^ The 
CSM has been used to correlate deviations from ideal 
tetrahedra in silicates as a function of applied pressure,^ 
and to analyze second-order Jahn- Teller systems across 
a variety of crystal structures.^ In these cases, the CSM 
was considered for the average crystallographic stucture, 
e.g. one where polyhedra possess a single CSM value. 
Here, we extend CSM to large-box modeling by calculat- 
ing it for every 0'Bi4 tetrahedron in the RMC supercell, 
obtaining distributions, rather than single values. 

The key finding to emerge from this study is that dis- 
placements from ideal atomic positions in Bi 2 Ti 2 OeO' 
and in particular, the nature of the 0'Bi4 tetrahedra in- 
dicate a tendency for Bi to lie in a disordered ring around 




FIG. 1: (Color online) The Bi 2 Ti 2 6 0' crystal structure is 
shown in (a), with 50% thermal ellipsoids representing the 
atomic displacement parameters from Rietveld refinement at 
14 K. The two sublattices are corner-sharing TiOe octahedra 
(blue) and corner-sharing 0'Bi4 tetrahedra (orange). Bi and 
O' are on ideal positions. Bi cations (black) appear as discs 
due to their displacement normal to the O'-Bi-O' bond. This 
disorder can be modeled using six-fold-split Bi and four-fold- 
split O' as shown in (b) (only the 0'Bi4 is sublattice shown). 



the ideal position, with some preference for near-neighbor 
Bi-Bi ordering. This reaffirms the case that, even when 
probed microscopically, Bi 2 Ti 2 OeO' is ice-like in its dis- 
order. We emphasize that in drawing the analogy with 
ice, we do not suggest the existence of ice-rules of the 
Bernal-Fowlei^ type in these systems. 



II. METHODS 

Synthesis and a detailed average structural analysis of 
the sample used in this study (including verification of 
purity) has been reported by Hector and Wiggin. 9 Briefly, 
a basic solution of titanium metal with hydrogen peroxide 
and ammonia was added to an acidic solution of bismuth 
nitrate pentahydrate and nitric acid. The resulting pre- 
cipitate was filtered, washed with a dilute ammonia solu- 
tion, dried at 50°C, and calcined in air for 16 h at 470°C. 
Time-of-flight (TOF) neutron powder diffraction on sam- 
ples held in vanadium cans was collected at the NPDF 
instrument at Los Alamos National Laboratory at 298 K 



and 14 K. NPDF is designed to collect high-resolution, 
high-momentum-transfer data suitable for production of 
the PDF, as well as traditional Rietveld refinement. We 
performed Rietveld refinement using GSAS.P Extrac- 
tion of the PDF with PDFGetI\P used Q max = 35 A -1 , 
and least-squares refinements of the PDF were performed 
with PDFguiP 

First-principles density functional methods were used 
to identify a possible ordered ground state. The lo- 
cal stability of Fd3m Bi 2 Ti 2 0y was investigated us- 
ing projector augmented wave potentials within the lo- 
cal density approximation as implemented in the vasp 
prograrrpmm and described previously.^ We found three 
zone-center lattice instabilities-one ferroelectric and two 
antiferroelectric-like modes-which were used to guide a 
systematic search for low-symmetry phases. We per- 
formed a series of structural relaxations with in each of 
the highest-symmetry isotropy subgroup d 39 * 40 ! consistent 
with the freezing-in of one or more of these lattice in- 
stabilities. We found that the lowest energy structure 
is ferroelectric (with a substantial polarization of P « 
20/iCcm~ 2 ), in the polar monoclinic space group Cm, 
consistent with the freezing-in of all three modes. 

Maximum entropy method (MEM) calculations of 
the nuclear scattering density were performed using 
PRIMAP RMC simulations were performed using 
RMCProfikP version 6ona5x5x5 supercell (cubic, 
sa 52 A per side) with 11000 atoms. These simulations 
were constrained by the PDF (in the form of £)(r))p3 
up to r — 12 A and by the Bragg profile of NPDF bank 
1, which contains the lowest-Q Bragg reflections. Crys- 
tal structures are visualized using VESTAP and Atom- 
Eye.^ Quantitative analyses shown in this work are 
taken from averaging many simulations in order to ob- 
tain an unbiased interpretation of the fit to data. Hard- 
sphere cutoffs were employed in RMC simulations to en- 
sure that atoms did not approach closer than the speci- 
fied nearest-neighbor distances, but no bunching was ob- 
served at these cutoffs. CSM for 0'Bi4 tetrahedra were 
calculated using a distance measure program provided by 
M. Pinsky and D. Avnir. 



III. RESULTS AND DISCUSSION 

A. Average structure from Rietveld refinement 

The profile resulting from TOF neutron Rietveld re- 
finement is displayed (for a single bank of data) in Fig. [2] 
The goodness of fit parameter R wp — 2.9%. The sample 
is the same as used in the Rietveld study of Hector and 
Wigging and the fit converges to similar structural pa- 
rameters (within error), so these will not be reproduced 
here. 

In Bi 2 Ti 2 OgO', preliminary Rietveld refinement using 
an ideal pyrochlore model can obtain a good fit to data, 
but large atomic displacement parameters (ADPs) indi- 
cate that significant atomic disorder is present, character- 




FIG. 2: (Color online) Neutron TOF Rietveld refinement 
of Bi2Ti20eO' at 14K (a) using the model of Hector and 
WiggirP gives a fit with R wp = 2.9%. Small features in the 
difference curve (b) at low Q are caused by shoulders of dif- 
fuse intensity around Bragg peaks. This component of the 
data can be interpreted using total scattering analysis. 



ized by an elliptical spread around the ideal atomic posi- 
tions. The ideal Bi position on the 16c (0,0,0) site can be 
fit using a large Uo 3 component. These appear as large, 
flat discs in Fig.Fna) with a radius of about 0.4 A. The 
discs envelop the Bi displacive disorder which we seek to 
accurately describe. The O' atoms at the 8a (§,§,§) 
position have large, isotropic displacements as well, cor- 
responding to Ui SO w 0.4A 2 . The TiOg sublattice, on the 
other hand, is described by small Ui SO values and does 
not display any signs of displacive disorder. Here we in- 
vestigate the precise nature of atomic displacements in 
the 0'Bi 4 sublattice. 

Improved refinement of the average pyrochlore struc- 
ture has been achieved by using a split-atom m odel 
for the A sites, such as in studies of Bi 2 Ti 2 060')2^ 
Bi 2 Sn 2 07Por La 2 Zr 2 07.^ Hector and Wiggin modeled 
Bi using a six-fold ring in the 96g position, 9 but acknowl- 
edged that their refinement does not clearly show a pref- 
erence for 96g versus 96/i (rotated 30° to each other) . The 
96g split-atom configuration is shown in Fig.[ljb). Com- 
parison to an ideal-position model in Fig.JlJa) shows that 
the Bi split-atom sites lie inside the anisotropic discs. 
Most O' intensity is still centered on the 8a site, but 
some occupancy is shifted away in 4 directions to form a 
tctrapod, modeled by partial occupancy of the 32e sites. 
From Rietveld refinement, we find that the Bi displace- 
ment parameters at 14 K and 300 K are similar in orien- 
tation and magnitude. This suggests that Bi displace- 
ments are frozen at room (and higher) temperatures, 
and what is being monitored in the scattering is not a 



snapshot of dynamic motion, but rather a description of 
static positions. The suggestion of frozen displacements 
at room temperature is consistent with measurements of 
the temperature- and field-dependence of the dielectric 
constant in Bi 2 Ti 2 6 0' thin films.^ 

Kinks in the difference profiles of the Rictvcld refine- 
ment in Fig.[2jb) can be attributed to shoulders of diffuse 
intensity around low-Q peaks. The diffuse scattering can- 
not be interpreted here because the Rietveld technique 
only models Bragg intensity. Crucial approximations are 
made to model a structure using only Bragg peaks: we 
must average any atomic correlations or distortions that 
do not possess long-range ordering. The Bi distortions 
are incoherent, and their description will require an ex- 
amination of diffuse scattering. Total scattering analysis 
of the real-space PDF, discussed here subsequently, pro- 
vides a real-space tool to model both Bragg and diffuse 
scattering simultaneously. 



B. Maximum entropy method 

The maximum entropy method (MEM) is used here to 
produce a map of neutron scattering density in the unit 
cell. The method, proposed originally by JaynesJ^ESl 
uses some testable information (in this case, observed 
Bragg peak intensities in TOF neutron diffraction), and 
is based on the most probable distribution (nuclear scat- 
tering density) being the one with the largest inform ation 
entropy. This method is described by SakatcP^I anc | i n 
the context of the MEM software program MEED,^ a 
precursor to the software code PRIMA used here. The 
MEM input consists of the observed Bragg structure fac- 
tor F b s (from a Le Bail fit, with phase information from 
the ideal structure) for a list of hkl reflections, the unit 
cell dimensions, space group, and the sum of all scat- 
tering lengths in the cell. The MEM as employed here 
is largely model-free in the sense that atomic positions 
are not specified during the calculation, but phase in- 
formation for F b s is biased by the ideal structure. The 
final result is required to obey the symmetry of the space 
group. This is also an average structure probe-diffuse 
scattering intensity is ignored. Thus all Bi atomic posi- 
tions in Bi2Ti20gO', for example, will be equivalent to 
each other as required by Fd3m. 

Despite these constrains, the MEM affords an excellent 
view of average Bi displacements. No prior description of 
Bi displacements, or even a knowledge of their existence, 
is used to produce them in the isosurface neutron den- 
sity unit cell displayed in Fig.[3^a). Arrows point to the 
circular Bi density that forms a ring around the O'-Bi- 
O' bond. Two-dimensional slices of the cell viewed along 
and normal to the O'-Bi-O' bond in Figs.[3^b) and[3|c), 
respectively, show additional detail of the Bi nuclear den- 
sity. The Bi ring appears hexagonal in Fig.[3^b), which 
is the required symmetry of the position in Fdim. The 
strongest intensity of Bi points toward the in-plane Ti 
in Fig.[3]jb). This corresponds to a 96/i site for Bi, not 
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FIG. 3: (Color online) The MEM nuclear density map of the 
Bi2Ti20eO' unit cell at 14 K is plotted in (a). Arrows point 
to Bi positions where the rings are rotated with the O'-Bi-O' 
direction almost parallel to the viewing direction. The circu- 
lar Bi isosurface is magnified in the inset. The 2-D density 
slice in (b) is the (111) plane, normal to O'-Bi-O' bonds, 
and (c) is the (121) plane with the O'-Bi-O' bond dashed. 
The imposed sixfold symmetry of the Bi ring is evident in the 
(111) slice. 

96<7p^l While the Bi shapes found by MEM agree with 
those from Rietveld refinements, there is no evidence of 
a tetrapod-shaped spread in the O' density in MEM. This 
shape would be allowed by symmetry. 

C. Least-squares PDF refinement 

When displacements of atomic positions and their cor- 
relations have been proposed in the pyrochlore struc- 
ture, it has been compared to the Si04 tetrahedral 
network of /3-cri stobalit e, which is isostructural to the 
0'Bi4 sublatticeP^MEEl This configuration would have 
long-range zig-zag ordering of Bi displacements shown in 
Fig. [5] These distortions may be ordered along any of 
the three Bi-Bi directions in the lattice, and the distor- 
tions in o ne dir ection need not be correlated with those 
in another t 23 ! 57 ! We use least-squares PDF refinement (a 
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FIG. 4: (Color online) Least-squares refinements of the 
Bi 2 Ti 2 060' PDF at 14 K. Panel (a) shows the fit to an ideal 
average structure (Bi centered on 16c sites with disc-shaped 
thermal parameters) as shown in Fig.JTJa). Panel (b) is the 
Cm structure from first-principles . The \ 2 is shown for two 
fits per model: a fit up to 3.4 A and a fit from 3.4 A to 6 A. 
Difference curves are shown below. 
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FIG. 5: (Color online) In (a), the ideal 0'Bi4 sublattice (or- 
ange tetrahedra) with split-atom Bi positions (black) can co- 
operatively order in a zig-zag (arrowed), /3-cristobalite-like 
pattern ordering seen as solid black lines. We use a double- 
tetrahedral model viewed along the O'-O' direction in (b) to 
investigate local correlations. This shape is viewed in (c) near 
the Bi-Bi direction with the most-probable wedges of the cir- 
cular Bi ADPs shown schematically for zig-zag ordering. 



D. Reverse Monte Carlo simulations 



small-box technique) to compare the observed local struc- 
ture to two models: the Fd3m Rietveld-derived average 
model with large ADPs, and the Cm ordered structure 
from first-principles calculations. 

The average model fit to the experimental PDF is 
shown in Fig.Qa). The large, disc-shaped ADPs on the 
Bi positions cause the peak at 2.2 A to flatten and dis- 
appear, and the fit does not significantly improve as r 
increases. The short-range fit up to r = 3.4 A and the 
medium-range fit from 3.4 < r < 6 A both give high 
X 2 values. The polar, ordered model fit in Fig.j4|b) de- 
scribes the low-r region very well. Here, a single Bi-O' 
distance is well-defined, so the peak at 2.2 A appears. 
The medium-range fit is poorer than the short-range fit, 
but still much better than the average model. From the 
Rietveld, MEM, and least-squares PDF analyses we know 
that the structure of Bi 2 Ti 2 060' shares some attributes 
of the average structure (lattice parameters, averaged po- 
sitions) and the local Cm model (bond distances and an- 
gles at low r). Analysis by least-squares refinement of the 
PDF is limited because we must construct a model that 
resembles both Fd3m and Cm. Due to the complex dis- 
order present in Bi2Ti2 0eO', we conduct reverse Monte 
Carlo simulations, where large-box models are generated 
by fits to the data (starting from the ideal structure), 
and extract statistical measures of these models. 



RMC simulations utilize large-box modeling with pe- 
riodic bo undary conditions and are not constrained by 
symmetry! 21 ! 42 ! 5 ^ This provides two principal advantages 
when simulating crystalline materials: the ability to 
model nuclear positions with arbitrary shapes, and the 
ability to investigate correlations between atoms on the 
A length-scale. We profit from both when modeling 
Bi2Ti2 0gO'. First, the true shapes of Bi and O' dis- 
placements must be elucidated. Second, Bi-Bi nearest- 
neighbor correlations may lead to cooperative distortions 
of the 0'Bi4 sublattice. This type of short-range order 
would be a signature of local regions where the lone-pair- 
active Bi displacements effect polar domains. 



1. Average atomic displacements 

The real-space local structure (coordination environ- 
ments, short-range correlations, etc.) of the RMC super- 
cell are driven by the fit to the PDF, while the Bragg 
fit constrains the long-range periodicity of the structure 
(or lack thereof) and ensures reliable displacement pa- 
rameters. In the case of Bi 2 Ti20eO', we monitor the 
agreement of Rietveld and RMC ADPs (calculated as the 
mean-square displacement from ideal positions). RMC 
simulations were started with the ideal pyrochlore lattice 
with all Bi atoms on 16c positions. As the simulation 
progresses, Bi atoms are always observed to move off the 
central position to form a ring that resembles the Ri- 



6 




^ 20 

in 

I 10 



I 1 I 1 

" (b) 


I 1 I 1 I 

• data 
— RMC 







Q (A j 

FIG. 6: (Color online) RMC fits to 14 K Bi 2 Ti 2 B 0' experi- 
mental data shown for (a) the PDF G(r) and (b) the Bragg 
profile. The RMC simulations employ supercells (5x5x5, 
« 52 A per side) that fit both data sets simultaneously. 
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etveld split-atom model. The fits to the 14 K PDF and 
Bragg profile are shown in Fig.[6]ja,b). 

A comparison of the ideal unit cell with a folded RMC 
supercell is shown in Fig. [7] In Fig.JT^b) we fold the 125 
unit cells of a supercell into a single unit cell. The result 
is a cell with "point clouds" on each atomic position that 
represent a map of the nuclear scattering density on each 
site in Bi 2 Ti 2 060'. Of particular note is the large spread 
of Bi and O' point clouds in Fig.[7]jb) in comparison to 
those of Ti and O. These distributions agree quantita- 
tively with the ADPs from Rietveld refinement. 

The 0'Bi4 sublattice from a 14 K RMC simulation is 
shown in Fig.^c), viewed close to the (110) direction. 
From this viewpoint, Bi rings are evident. These rings 
open normal to the O'-Bi-O' bonds, as expected from 
Rietveld ADPs and MEM. The average shape of a Bi ring 
is seen in Fig.[8]ja). Viewed along the O'-Bi-O' direction, 
it appears circular with a radius of about 0.4 A. A his- 
togram of Bi distances from the ideal 16c site is given in 
Fig.[8^c), with a clear maximum at the ring radius. The 
same procedure applied to O' yields a tetrapodal object. 
Seen in Fig.^b), the arms of this tetrapod point away 
from the four neighboring Bi atoms. One arm of the 
tetrapod is in the center of the plot, pointing normal to 
the page. A key distinction between O' and Bi is that 
the dense cluster of O' positions is still centered on the 
ideal 8a position, seen as a bright cluster in the center of 
Fig.|8^b). In contrast, Bi intensity is low at the center and 
most intensity lies on the ring perimeter. The displace- 
ment histograms in Figs.[8|c,d) reveal similar shapes in 
RMC simulations of 14 K and 300 K data, with a slight 
broadening at high temperature. This supports the idea 



FIG. 7: (Color online) The ideal, isotropic Bi 2 Ti20 6 0' unit 
cell is plotted in (a) along the a direction, with Bi atoms in 
black, Ti blue, and O orange. This unit cell can be visually 
compared with the folded 14 K RMC supercell in (b). Posi- 
tional disorder is much greater in Bi and O' (which appear 
X-shaped in this projection) than in Ti and O (appear as 
clusters of four O). In (c), the 0'Bi4 sublattice is plotted near 
(110) to show the Bi rings. 



that these displacements are frozen at high temperatures 
and do not fundamentally change with cooling. 

The RMC model corroborates with Rietveld and MEM 
while providing further evidence for a tendency toward a 
sixfold Bi ring. In Fig.^a), a histogram of the the rota- 
tion angle 9 of Bi around the (111) axis is fit to a sixfold 
cosine curve. This curve, with a minimum at 9 = 0°, 
corresponds to the 96g sites and does not fit the data. 
Instead, Fig.^b) shows excellent agreement between the 
RMC result and the 96/i six-fold curve, which is shifted 
by 30° from 96<?. Note however that the intensity at min- 
ima in the curves is not zero; there is still some tendency 
for Bi to lie at any 9. It is not known whether increased 
preference for the 96/i site would push the minima to be 
nearly zero. It is therefore possible that stronger ordering 
exists, but is obscured by the resolution of the data. 



2. Correlated distortions: double-tetrahedral model 

Cooperative /3-cristobalite correlations can be visual- 
ized in the polyhedral configuration given in Fig.[5j The 
motif of two corner-linked 0'Bi4 tetrahedra is the basis 
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FIG. 8: (Color online) Atomic positions for (a) Bi and (b) O' 
at 14 K from 14 independent RMC simulations are plotted in 
the {111} plane normal to the O'-Bi-O' bond. Histograms 
of displacement distances are plotted for 14 K and 300 K 
simulations in (c) and (d), with the 14 K histograms shaded 
for clarity. The Bi median displacement of 0.4 A agrees with 
the Bi Rietveld ADP. 
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FIG. 9: Histogram of the average Bi displacement angle 8 
around the O'-Bi-O' bond at 14 K. Solid lines show the ex- 
pected six-fold cosine curves that indicate preference for (a) 
96g or (b) 96/i hexagonal ring symmetry. 



of our analysis. The zig-zag pattern is comprised of Bi 
that are all in the same plane (denoted with arrows). 
This plane is defined by the displacement of any Bi atom 
in the chain, so the central Bi atom can only partici- 
pate in zig-zag ordering in one direction at a time. The 
two inline nearest neighbors of the central atom should 
have a tendency to displace opposite from the central 
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FIG. 10: Schematic of possible values of <f>, probing nearest- 
neighbor correlations of Bi displacements (not to scale). 
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FIG. 11: (Color online) Histograms (a) of the Bi-Bi angle cj> 
for neighboring displacements, with error bars of ±<r. The dif- 
ference between the RMC result and the same configurations 
with randomly shuffled Bi displacements is shown in (b). The 
difference curve peaks at (f> = 180°, indicating some zig-zag- 
like correlation. Angles around <f> = 120° seem unfavorable. 



Bi displacement vector. We can therefore quantify the 
correlation by examining the angle <j> between the dis- 
placement vector of the central Bi versus the displace- 
ment vectors of the two relevant nearest neighbors. Each 
double-tetrahedral shape in the supercell is examined in 
this manner. A schematic showing selected orientations 
of neighboring Bi displacements and their corresponding 
values of </> is given in Fig.flO) 

In Fig.[llja), the relative angle <j> histogram is plotted 
for the RMC simulations. It does not have the same six- 
fold modulation as the averag ed Bi angles 6 in Fig.^b) 
because cj> is defined relative to its neighbors, while 9 is 
defined relative to the crystal axes. We investigate the 
tendency for ordering by shuffling Bi displacements: the 
set of all individual Bi displacements from their ideal 
sites is preserved, but redistributed randomly among the 
Bi atoms. As a result, their orientations with respect 
to each other are disrupted. Note that, as a result of 
the removal of this local correlation, the shuffled <j> curve 
is simply a six-fold cosine curve with modulation analo- 
gous to 6 in Fig.[9jb). The difference between the RMC 
and shuffled <f>, plotted in Fig.[TT|b), shows a peak at 
<f> = 180°, indicating a preference for zig-zag ordering. 
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FIG. 12: Histogram (a) of the CSM tetrahedricity S(T d ) for 
three supercells: the RMC supercell, the same RMC supercell 
with Bi displacements randomly swapped, and a model with 
random circular Bi rings of equal radius. The peak at higher 
S(Td) in the circular model indicates more deviation from 
ideality than in the RMC supercells. These histograms are 
plotted as cumulative probability functions in (b), where the 
agreement between the RMC and swapped supercells suggest 
little evidence for cooperative short-range Bi ordering within 
tetrahedra. A fully-ordered model is denoted by a vertical 
line at S{T d ) = 0.008. 



Correlations where 60° < <j> < 120° are not preferred. 
While these results point to anti-alignment of Bi displace- 
ments, Bi 2 Ti 2 060' should not be considered to follow 
the /3-cristobalite model strictly, since Bi would need to 
prefer the 96g site rather than 96h. 



line in Fig. 12 a), and again as a cumulative fraction in 
FigQb). 

In a case where zig-zag (or otherwise correlated) dis- 
tortions are present on most 0'Bi4 tetrahedra in the 
supercell, we expect Bi displacements to be correlated 
with their nearest neighbors, and a polyhedral rigidity 
reminiscent of Si04 tetrahedra might exist. If these dis- 
placements from ideal positions were randomly swapped 
among Bi atoms, the ordering would be destroyed and 
S(T ( i) would increase toward the random circular model. 
In a case where there is no short-range Bi-Bi ordering, 
randomly swapping the displacements should have no ef- 
fect on the S(Td) distribution. 

In Fig. [l2] we compare the S(T d ) distributions for the 
RMC supercell and the same supercell with Bi displace- 
ments randomly swapped. The distributions are shown 
as histograms of S(T d ) in (a) and cumulative fractions 
in (b). The random swapping slightly increases S(T d ) 
for the RMC supercell, seen as a shift to the right in 
Fig.[l2|b). This curve is still quite separated from the 
random circular model. 

Two aspects of the Bi distribution can be gleaned from 
this comparison. First, the Bi displacements of the RMC 
supercell are more strongly tetrahedral than the random 
circular model. This could be due to rigidity of the 0'Bi4 
tetrahedra, or the fact that the RMC has a distribution 
of Bi displacements [Fig.[8^c)] whereas the circular model 
has strictly r — 0.4 A. The separation between RMC and 
shuffled CSM points to rigidity that may be explained by 
correlated Bi displacements. Just as shuffling Bi displace- 
ments disrupts the ordering between pairs of Bi as viewed 
by (j) in Fig.[TT|(b) , it also disrupts ordering between sets 
of four Bi displacements as viewed by CSM in Fig. 12 'b). 



IV. CONCLUSIONS 



3. Correlated distortions: continuous symmetry measures 

The CSM provides a quantitative measure of a poly- 
hedron's distortion, in the form of a "distance" from 
idealityPfiKsH j n our case, a perfect tetrahedron would 
have tetrahedricity S(T d ) = 0, while distortions from 
ideality increase this value toward 1. As a point of ref- 
erence, the Jahn- Teller distorted CUO4 tetrahedron in 
CuMn 2 0j2^ has S(T d ) = 0.0235 while a square plane 
has S(T d ) = §. 

Application of CSM to the 0'Bi4 sublattice is useful 
because S(T d ) describes the correlations between four 
neighboring Bi atoms. Comparison between models and 
RMC simulations are shown in Fig. [12] For a model with 
ordered /3-cristobalite-type ordering with a similar dis- 
placement magnitude, S(Td) = 0.008. This is plotted 
as a vertical line in Fig. 12jb). In contrast, uncorrelated 
and evenly circular (with radius 0.4 A) Bi displacements 
produce a wide distribution of S(T d ), shown as a dotted 



We have shown that a structural investigation of the 
strongly disordered pyrochlore Bi 2 Ti 2 OeO' benefits from 
modeling real- and reciprocal space structure simulta- 
neously. The combination of Rietveld refinement, MEM, 
least-squares PDF refinement, and RMC simulations pro- 
vide a description of the structure over many length 
scales. Of particular interest is the behavior of the dis- 
placed lone-pair active Bi 3+ cation. Average structure 
models indicate a preference for static displacements. 
MEM and RMC show that Bi prefers the 96h (point- 
ing between nearby O from the TiOe sublattice) rather 
than 96g sites (pointing toward neighboring O). Statis- 
tical measures of these displacements are extracted from 
RMC supercells. 

RMC simulations show that Bi displacements lie at all 
angles 9 in a ring normal to the O'-Bi-O' bond with ra- 
dius r ~ 0.4 A, but there is a preference for 9 correspond- 
ing to 96h positions. O' atoms are centered on the ideal 
site, but in a tetrapodal shape with arms pointing away 
from neighboring Bi. Comparing the RMC result with an 
identical supercell with shuffled Bi displacements reveals 
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that neighboring Bi displacements show some correlation 
with each other, with a preference for phi = 180° align- 
ment. This corresponds to a zig-zag type ordering that 
merits further investigation. Measures of the tetrahedric- 
ity S{Td) of the 0'Bi4 sublattice indicate that neighbor- 
ing displacements are weakly correlated, thus the sublat- 
tice may exhibit some rigidity. 
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